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• 1—1 ■ Abstract 

Q_^| of perfect, or ideal, conductors. It is analogous to Thomson's theorem 

on the equilibrium electric field and charge distribution in a system 

of conductors. We first prove Thomson's theorem using a variational 

^ I principle. Our new theorem is then derived by similar methods. We 

find that magnetic energy is minimized when the current distribution 
If^ . is a surface current density with zero interior magnetic field; perfect 

'^ I conductors are perfectly diamagnetic. The results agree with currents 

•^ ' in superconductors being confined near the surface. The theorem im- 

plies a generalized force that expels current and magnetic field from 
the interior of a conductor that loses its resistivity. Examples of solu- 
tions that obey the theorem are presented. 
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1 Introduction 

Thomson's theorem states that electric charge density on a set of fixed con- 
ductors at static equihbrium is distributed on the surface of the conductors 
in such a way that the interior electric field is zero and the surface electric 
field is normal to the surface. Here we will prove an analogous theorem on 
the magnetic field and current distribution in ideal conductors. We find that 
a stationary current density must distribute itself on the surface of the con- 
ductors in such a way that the interior magnetic field is zero while the surface 
magnetic field is perpendicular to both the current density and the surface 
normal. 

When W. Thomson (Lord Kelvin) derived his theorem in 1848 a mag- 
netic analog did not seem interesting since conductors with zero resistivity 
were unknown. Since the discovery of superconductivity in 1911 this has 
changed and stationary current distributions in type I superconductors, be- 
low the critical field, indeed obey our theorem. In spite of this the theorem 
is not stated, hardly even hinted at, in the literature. We wish to emphasize 
that even though the theorem applies to superconductors, phase transitions, 
statistical mechanics, or thermodynamics are irrelevant. Zero resistivity is 
simply assumed, not explained or derived. The theorem is purely a conse- 
quence of classical electromagnetism. 

The outline of this article is as follows. We start by deriving Thomson's 
theorem using a variational principle. We then derive our minimum mag- 
netic energy theorem in an analogous way and discuss previous work on the 
problem. After that an illuminating example is presented in which the en- 
ergy reduction due the interior field expulsion can be calculated explicitly. 
An Appendix gives further motivation and explicit solutions illustrating the 
theorem for simple systems. 

2 Energy minimum theorems 

Electromagnetic energy can be written in a number of different ways. Here 
we will assume that there are no microscopic dipoles so that distinguishing 
between the Z), H and E, B fields is unnecessary. That this is valid when 
treating the Meissner effect in type I superconductors is stressed by Carr [1]. 



Relevant energy expressions are then, 

E, + E^ = ^f{E^ + B^)dV (1) 

on J 
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Here the first form is always valid while the two following assume quasi 
statics, i.e. essentially negligible radiation. 

According to Thomson's theorem, the electric charge on a set of con- 
ductors distributes itself on the conductor surfaces thereby minimizing the 
electrostatic energy. W. Thomson did not present a formal mathematical 
proof but such proofs may be found in most classic textbooks [21 [3l HJ |5] . A 
recent derivation of the theorem in its differential form is by Bakhoum [B]. 
A derivation based on a variational principle can be found in the textbook 
by Kovetz [7]. A different approach also based on a variational principle, is 
presented below. Thomson's result is widely known and is useful in many 
applications. It has e.g. been used to determine the induced surface charge 
density [SI E] , and in the tracing and the visualization of curvilinear squares 
field maps [nj- Other applications range from interesting teaching tools [10] 
to useful computational methods such as Monte Carlo energy minimization 

There are various similarities between electrostatics and magnetostatics, 
or quasi- statics, but for resistive media the magnetic field due to current dis- 
sipate^j. For perfect., or ideal conductors, however, there should be something 
corresponding to the magnetic version of Thomson's theorem. Indeed, below 
we will prove a theorem analogous to that of Thomson: Magnetic energy 
is minimized by surface current distributions such that the magnetic field is 
zero inside while the surface field is normal to the current and the surface 
normal. Energy conservation is assumed restricting the validity to perfect, 
or ideal, conductors. Previously somewhat similar results have appeared in 
the literature [121 US] and we discuss those below. 



^Note that we are not concerned with magnetism due to microscopic dipole density. 



2.1 Thomson's theorem 



In equilibrium, the electrostatic energy functional for a system of conductors 
surrounded by vacuum, may be written acl. 
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by combining the electric parts of ([I]) and (^ and using E = —Vcj). We 
now split the integration region into the volume of the conductors, V^n, the 
exterior volume, Kut, and the boundary surfaces S, 
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where a is the surface charge distribution. 

We now use, (V0)^ = V ■ (0V0) — 0V^0, and rewrite the divergencies 
using Gauss theorem. The energy functional then becomes. 
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where V"*" and V~ are the gradient operators at the surface in the outer and 
inner limits, respectively. The total charge in each conductor is constant and 
restricted to the conductor volume and surface. We handle this constraint 
by introducing a Lagrange multiplier A. Infinitesimal variation of the energy 
then gives. 
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From this energy minimization, the Euler-Lagrange equations become: 
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^Thc ultimate motivation for this specific form of tlie energy functional, and the cor- 
responding one in the magnetic case, is that they lead to simple final equations. 
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According to equations (|H]) the potential is constant inside the conductor 
in the niinimum energy state, and therefore ^ = there. Equations ([H]) 
mean the electric charge is distributed on the surface in such a way that the 
potential is constant there. The second of eqs. (|8]) implies that V^0 = and 
the first of eqs. (19]) then implies that n ■ E^ = Ana. This concludes the proof 
of Thomson's theorem. 



2.2 Minimuin magnetic energy theorem 

A similar procedure will now be applied to the magnetic field. We write the 
magnetic energy functional for a time independent magnetic field as, 
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i.e. as two times the form ([2]) of the magnetic energy minus the form ([T]), 
using B = V X A. As before we split the volume into the volume interior 
to conductors, the exterior vacuum, and the surface at the interfaces, and 
write, 
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where k is the surface current density. We now use the identity, 
(V X A)^ = V ■ [A X (V X A)] + A ■ [V X (V X A)] , 
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and then use Gauss theorem to rewrite the divergence terms. The energy 
functional then becomes: 
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As in the electric case, constraints must be imposed. Due to charge conser- 
vation, the electric current density must obey the continuity equation, for 
zero charge density, both inside and on the surface of the conductors 
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where Vs is the surface gradient operator. Notice that these constraints are 
local, not global. In other words, the relevant Lagrange multiplier is not 
constant but a scalar field A(r). Using this, infinitesimal variation of the 
magnetic energy gives. 
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Equating this to zero we find that, 
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Kut : V X (V X A) = V X S 

are the Euler-Lagrange equations for this energy functional. 

According to eq. flTSj) B = V x VA = 0, so the magnetic field must be 
zero in Vin- Consequently also the volume current density is zero, j = 0, 
inside the conductor, in the minimum energy state. 



2.2.1 Surface currents 

Now consider the results for the surface, eq. (TT^ . Our results from V^n show 
that V^ X A = 0, so the equation reads. 
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V B^O j = 

vacuum 




y ideal conductor 
'° B^O j = 

Figure 1: Some of the results of our theorem on the current density and magnetic 
field of an ideal conductor at minimum energy are illustrated here. Here B^ is 
the magnetic field on the outside of the surface S with surface unit normal n. The 
bulk current density j is zero, only the surface current density k is non-zero. 



Let us introduce a local Cartesian coordinate system with origin on the sur- 
face, such that the surface is spanned by x, y with unit normal fi = z = xxy. 
Assuming that the surface is approximately flat we then have that. 
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Since, A = V^A, the vector potential is tangent to the conducting surface 
and we get. 
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for the surface current density. Rewriting the triple vector product we find. 
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so the surface current density is parallel to the outside normal derivative of 
the vector potential. We note that this agrees with the well known result [S], 
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for the case of zero interior field {B_ = 0). Our results are summarized in 

Fig.m 

2.2.2 External fields 

Neither Thomson's theorem nor our minimum magnetic energy theorem are 
formally valid for conductors in constant external fields. In both cases, how- 
ever, such a situation can be regarded as a limiting case. In the case of Thom- 
son's theorem one can include two large, distant, and oppositely charged par- 
allel conducting plates. A small system of conductors between these can then 
approximately be regarded as in an electric field that approaches a constant 
external field at large distance. In a similar way the set of perfect conductors 
can be thought of as inside two large perfectly conducting Helmholtz coils 
(tori) which provide an approximately constant external magnetic field at 
large distance. 

2.3 Previous work 

The fact that there is an energy minimum theorem for the magnetic energy of 
ideal, or perfect, conductors, analogous to Thomson's theorem, is not entirely 
new. In an interesting, but difficult and ignored, article by Karlsson [12] such 
a theorem is stated. Karlsson, however, restricts his theorem to conductors 
with holes in them. In the electrostatic case charge conservation prevents the 
energy minimum from being the trivial zero field solution. In our magnetic 
ideal conductor case the corresponding conservation law is the conservation 
of magnetic fiux through a hole fi5\ [T6] . As long as one conductor of the 
system has a hole with conserved fiux there will be a non-trivial magnetic 
field. To require that all conductors of the system have holes, as Karlsson 
does, therefore seems unnecessarily restrictive. One of Karlsson's results is 
that a the current distribution on a superconducting torus minimizes the 
magnetic energy. 

A result by Badia-Majos [T3] comes even closer to our own and we outline 
it here. The current density is assumed to be of the form, 

j = qnv, (27) 

where q is the charge of the charge carriers and n is their number density. 



The time derivative is then given by, 
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assuming that only the Lorentz force acts (ideal conductor). We now recall 
Poynting's theorem [17] for the time derivative of the field energy density of 
a system of charged particles, 
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The first term on the right hand side normally represents resistive energy 
loss. Here we use the result for E from eq. ( l28l) . 



E 



m dj 



1 



q'^n d t qnc 



JxB, 



and get. 



j-E 



m dj 



m 



■J 



(30) 



(31) 



q'^n dt dt \2q'^n J 

This is thus the natural form for this term for perfect conductors. We insert 
it into (129|) . neglect radiation, and assume that E'^ <^ B^. This gives us. 
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Finally inserting, j = (c/47r)V x B, here, gives. 
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for the conserved energy, after integration over space and time. 

Badia-Majos p3] then notes that this energy functional implies flux ex- 
pulsion from superconductors. Variation of the functional gives the London 
equation [T8] . 

B + ^^Vx{VxB) = 0. (34) 

Badia-Majos, chooses not to point out that this classical derivation of flux 
expulsion is in conflict with frequent text book statements to the effect that 
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no such classical result exists. Further work by Badia-Majos et al. [T^ on 
variational principles for electromagnetism in conducting materials should be 
noted. 

Finally we should mention the pioneering work by Woltjer [20j on energy 
extremizing properties of, so called, force free magnetic fields. In plasma 
physics a lot of further work has been done in that tradition. It has, however, 
not been concerned with currents and fields inside or on the boundary of 
bounded domains separated by vacuum, as we are here. 

3 Ideally conducting sphere in external field 

We now know that magnetic energy minimum occurs when current flows only 
on the surface and the magnetic field is zero inside. Let us consider a perfectly 
conducting sphere surrounded by a fixed constant external magnetic field. 
Here we will calculate the surface current needed to exclude the magnetic 
field from the interior and how much the total magnetic energy i^^ is then 
reduced. In order to exclude a constant external field B^. from its interior 
the currents on the sphere must obviously produce an interior magnetic field 
Bi = —Be, thereby making the total field B = Be + Bi zero in the interior. 
We will use that a constant field is produced inside a sphere by a current 
distribution due to rigid rotation of a constant surface charge density [21] . 



3.1 Energy of the external field 

For the total magnetic field to have a finite energy Em we can not assume that 
the constant external field extends to infinity. Instead of using Helmholtz 
coils to produce it we simplify the mathematics and produce our external 
field by a spherical shell of current that is equivalent to a rigidly rotating 
current distribution on the surface. This can be done in practice by having 
as set of rings representing closely spaced longitudes on a globe with the 
right amount of current maintained in each of them. Such a spherical shell 
of rigidly rotating charge produces a magnetic field that is constant inside 
the sphere and a pure dipole field outside the sphere (see Fig. |2]): 

——T- tor r < it 

Be(r) = I (35) 

3{m ■ f)f — m 






ioi r > R 
10 




Figure 2: The field lines of the field of eq. (I35p for a nonmagnetic sphere with a 
rigidly rotating homogeneous surface charge density. 

Here r = |r| and the center of the sphere is at the origin. If Q is the total 
rotating surface charge and a; its angular velocity, 
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see eq. f lM|) below. It is now easy to calculate the magnetic energy of this 
field. One findqfl, 
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Inside this sphere, which is assumed to maintain a constant current density 
on its surface, we now place a smaller perfectly conducting sphere. 



3.2 Magnetic energy of the two sphere system 

We assume that the small sphere in the middle of the big one has radius a < R 
and that it also produces a magnetic field by a rigidly rotating charged shell 
on its surface. We denote its dipole moment by rrii so that its total energy 
would be, 
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■^If ([36]) is inserted for m here we get i?°j = {uj/c)'^RQ'^/9 which is equal to £",„ of eq. 
56|) below, for ^ = 1, corresponding to surface current only, as it should. 
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if it was far from all other fields, according to our previous result f l37p . We 
now place the small sphere inside the large one and assume that rrii makes 
an angle a with m = mm, 

m ■ rrii = mmi cos a. (39) 

The total energy of the system is now, 

E^ = ^ /(Be + B.fdV = El + El^ + E„ (40) 

where the coupling (interaction) energy is, 

E, = ^JB,-B,dV. (41) 

This integral must be split into the three radial regions: 0<r<a, a<r< 
i?, and R < r. The calculations are elementary using spherical coordinates. 
The contribution from the inner region is, 

Eci = ^f B^.B,dV=^'^cosa. (42) 

4:71 Jr<a 3 R^ 

The middle region, where there is a superposition of a dipole field from the 
small sphere and a constant field from the big one, contributes zero: Ec2 = 0. 
The outer region gives Ec-i = [2 / 3)1711711 cos a /R^. Summing up one finds, 

^0 = 2^ cos a, (43) 

for the magnetic interaction energy of the two spheres. 

3.3 Minimizing the total magnetic energy 

The total magnetic energy of the system discussed above is thus, 

2 2 

in 171 TflTfl ' 

Era{m,,a) = El + El^ + E, = — + ^ + 2-^ cos a. (44) 

We assume that m and ttij are positive quantities. This means that as a 
function of a this quantity is guarantied to have its minimum when cos a = 
—1, i.e. for a = 7T. Thus, at minimum, the dipole of the inner sphere has the 
opposite direction to that of the constant external field, rrii = —rriirh. 
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Now assuming a = vr we can look for the minimum as a function of m,. 
Elementary algebra shows that this minimum is attained for, 

mi = i—j m = mimin- (45) 

The magnetic field in the interior of the inner sphere (r < a) is then, 

B, + B, = ^ + ^^^ = f ^ - ^^l rh = 0, (46) 

so it has been expelled. The minimized energy (jH]) of the system is found to 
be, 

31 
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The relative energy reduction is thus given by the volume ratio of the two 
spheres. 

Using E^ from fl37|) . the energy lowering is now found to be: 

This result is independent of the radius R of the big sphere introduced to 
produce the constant external field. It shows that the energy lowering corre- 
sponds to three times the external magnetic energy in the volume 47ra^/3 of 
the perfectly conducting interior sphere. 



4 The mechanism of flux expulsion 

In 1933 Meissner and Ochsenfeld |22], discussing their experimental discov- 
ery, stated that it is understandable that an external magnetic field does not 
penetrate a superconductocl but that the expulsion of a pre-existing field at 
the phase transition cannot be understood by classical physics. This state- 
ment has since been repeated many times. We are not aware, however, of any 
deeper investigations of what classical electromagnetism predicts regarding 
the behavior of perfect conductors in this respect, at least not prior to the 



^Eddy currents induced in accordance with Lenz law do not dissipate because of zero 
resistivity. 
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work of Karlsson [12] and Badia-Majos [13]. Our theorem strengthens the 
conclusion of Badia-Majos that a magnetic field is expelled according to clas- 
sical electromagnetism. Hirsch [23] has pointed out that the Meissner effect 
is not explained by BCS theory. Here we briefly speculate on the microscopic 
physical mechanism of this expulsion. 

Assume that a resistive metal sphere is penetrated by a constant magnetic 
field. Lower the temperature until the resistance vanishes. How does the 
metal sphere expel the magnetic field, or equivalently, how does it produce 
surface currents that screen the external field? According to Forrest [2l| this 
can not be understood from the point of view of classical electrodynamics 
since in a perfectly conducting medium the field lines must be frozen-in. This 
claim is motivated thus: When the resistivity is zero there can be no electric 
field according to Ohm's law, since this law then predicts infinite current. 
But if the electric field is zero the Maxwell equation, V x ^ — dB/dt = 0, 
requires that the time derivative of the magnetic field is zero. Hence it must 
be constant. 

This argument is flawed since Ohm's law is not applicable. Inertia, induc- 
tive or due to rest mass, prevents infinite acceleration. Instead the system of 
charged particles undergoes thermal fluctuations and these produce electric 
and magnetic fields. These fields accelerate charges according to the Lorentz 
force law. In the normal situation the corresponding currents and fields 
remain microscopic. When there is an external magnetic field present the 
overall energy is lowered if these microscopic currents correlate and grow to 
exclude the external field. According to standard statistical mechanics the 
system will then eventually relax to the energy minimum state consistent 
with constraints. We note that Alfven and Falthammar [25j state that "in 
low density plasmas the concept of frozen-in lines of force is questionable" . 

Another argument by Forrest [23] is that the magnetic flux through a per- 
fectly conducting current loop is conserved. Since one can imagine arbitrary 
current loops in the metal that just lost its resistivity with a magnetic field 
inside, the field must remain fixed, it seems. Is is indeed correct that the flux 
through an ideal current loop is conserved, but the actual physical current 
loop will not remain intact unless constrained by non-electromagnetic forces. 
There will be forces on a loop of current that encloses a magnetic flux that 
expands it [5]. This is the well know mechanism behind the rail gun, see e.g. 
Essen [26]. All the little current loops in the metal will thus expand until 
they come to the surface where the expansion stops. In this way the interior 
field is thinned out and current concentrates near the surface. So, the flux is 
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conserved through the loops, but the loops expand. 

It was noted early in the history of superconductivity that the Meissner 
flux expulsion is necessary if superconductors obey normal thermodynamics. 
Gorter and Casimir [27] observed that the flnal thermal equilibrium state 
of a superconductor must be independent of whether the external magnetic 
field existed inside the body prior to the phase transition or if it was added 
after the transition already had occurred. Our theorem indicates that it is 
perfectly natural that an interior magnetic field is expelled in the approach 
to thermodynamic equilibrium. 

5 Conclusions 

It is amazing that the theorem derived here has not been stated before. Zero 
resistivity conductors have been known since 1911 and zero resistivity is also 
considered a good approximation in many plasmas, so the ideal, or perfect 
conductor, is a well known concept and its magnetic energy minimum ought 
to be of great importance. It is clear from results above, and further moti- 
vated by the examples in the Appendix below, that type I superconductors 
below their critical field obey the theorem, and the reason that these only 
have surface current and zero interior field is thus simply minimization of 
magnetic energy. Naturally there is some other energy involved that is re- 
sponsible for the zero resistivity itself but apart from being implicitly assumed 
constant in our variations it is irrelevant to the current investigation. 

A Appendix 

In this Appendix we present further evidence in support of our theorem. The 
purpose of the calculations presented here is to illuminate and elucidate the 
physical meaning and the mechanisms behind the expulsion of current and 
magnetic field from the interior of ideal conductors. 

A.l Magnetic energy minimization in simple one de- 
gree of freedom model systems 

We investigate two simple one degree of freedom model systems and use them 
to illustrate how the minimum magnetic energy theorem works. We take 
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systems in which the magnetic energy Em can be calculated exactly so that 
energy minimization amounts to minimizing a function of a single variable. 
The systems are both related to a system used by Brito and Fiolhais p!0] to 
study electric energy. 

A. 1.1 Magnetic energy of coaxial cable 

The minimum magnetic energy theorem can be illustrated in such a simple 
system as a coaxial cable. The cable can be modeled by an outer cylindrical 
conducting shell with radius b, carrying an electric current /, and a concentric 
solid cylindrical conductor with radius a < b, carrying the same electric 
current in the opposite direction. We now assume that the total current on 
the inner cylinder is the sum of surface current Ig and bulk interior current 
4. See Fig. m 




Figure 3: Notation for the coaxial cable. Magnetic energy is minimized when the 
current on the inner conductor is pure surface current Is = I, and I^ = 0. 

Here we use cylindrical coordinates, p, ip, z and put Is = (1 — ri)I at p = a 
for the surface current, and I^ = rjl for the bulk current in < p < a. Then 
I = Is + Iv, and we get the magnetic field. 
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using Ampere's law. Thus, the magnetic field energy for a length L of the 
cable is given by: 
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This magnetic energy reaches its minimum for zero bulk current, t] = 0, 
corresponding to surface current only, and zero field for < p < a. 

A. 1.2 Current in sphere due to rigidly rotating charge 

Consider an ideally conducting sphere of radius R. Assume that there is a 
circulating current in the sphere which can be seen as the rigid rotation of a 
charge Q evenly distributed in the thick spherical shell between r = a < R 
and r = R. The charge density. 
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is assumed to rotate with angular velocity oj = u z relative to a an identical 
charge density of opposite sign at rest. The current density is then. 



3{r) = Q{r)u: x r. 



(52) 



and the current, / = ^w, passes through a half plane with the z- axis as 
edge. 

The vector potential produced by this current density can be found using 
the methods of Essen ^ , see also [221 EDI EI] • If we introduce ^ = a/R we 
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find, 
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for < r < a 



for a < r < i? (53) 



for R < r 



The parameter ^ = a/R is zero, ^ = 0, for a homogeneous ball of rotating 
charge, while ^ = 1 corresponds to a rotating shell of surface charge, see Fig. 
m Comparing with the vector potential for a constant field, A = ^Bq x r 
we see that, 

Quj {i-e) _Qu^ (1+0 ..,. 

° ci?(i-a ci?(i+e+a ^ ^ 

is the field in the central current free region < r < a. 




a=^R 



R 



Figure 4: Some notation for the system considered here. Current density flows 
in a thick spherical shell as a rigid rotation of constant charge density between 
r = a = ^R and r = R. 



A. 1.3 Magnetic energy of rotating spherical shell current 

We now calculate the magnetic energy of this system using the formula, 

E^ = ^Jj-AdV. (55) 

Performing the integration using spherical coordinates gives, 



i?r.(a;,0=(^)'^/(0, (56) 



where. 



f(^) - 35(1 + ^ + ^2)2 ' (57) 

is a function of the dimensionless parameter ^. Note that /(O) = 2/35, that 
/(I) = 1/9, and that /(^) is monotonically increasing, by a factor of almost 
2 in the interval to 1. 

This expression for the energy is the Lagrangian form of a kinetic energy 
which depends on the generalized velocity u = ip, 

Lm{v,0 = ^^f{Ov'- (58) 

Since the generalized coordinate if does not appear in the Lagrangian L^ 
the corresponding generalized momentum (the angular momentum), 

is a conserved quantity. The corresponding Hamiltonian, and relevant, ex- 
pression for the magnetic energy is then H^ = p^p'p — L^, expressed in terms 
of P^, 

EUO = HUp.,0 = Jq^^- (60) 

The function l//(0 is plotted in Fig. |5l We now consider the two energy 
expressions (1551) and fpU]) separately. 

Case of constant current: We first consider the case that the current, 
I = Quj/2tt, is constant. Changing ^ then means changing the conductor 
geometry while keeping a constant total current, or, equivalently, angular 
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Figure 5: A graph of the function l//(^) which is proportional to the Hamiltonian 
form of the magnetic energy (j60p of our model system. Note that ^ = corresponds 
to volume (bulk) current and ^ = 1 to pure surface current. 



velocity u = (p. One might regard the total current as flowing in a con- 
tinuum of circular wires. Changing ^ from zero to one means changing the 
distribution of these circular wires from a bulk distribution in the sphere to 
a pure surface distribution, while maintaining constant current. According 
to a result by Greiner [32] a system will tend to maximize its magnetic en- 
ergy when the conductor geometry changes while currents are kept constant. 
This has also been discussed in Essen |26]. In conclusion, if currents are kept 
constant the magnetic energy fl58l) will tend (thermodynamically) to a stable 
equilibrium with at a m,aximum, value and we note that this corresponds to 
a pure surface current ^ = 1. 

Case of constant angular momentum: Assume now that we pass 
to the Hamiltonian (canonical) formalism. Thermodynamically this type of 
system should tend to minimize its phase space energy (!60|) in accordance 
with ordinary Maxwell-Boltzmann statistical mechanics. As a function of ^ 
this Hamiltonian form of the energy Em{C) clearly has a minimum at ^ = 1, 
see Fig. O corresponding to pure surface current. In this case therefore there 
will be current density only on the surface in the energy minimizing state. 
This is in accordance with the our minimum magnetic energy theorem. It 
is notable that both the assumption of constant current and the assumption 
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of constant angular momentum lead to a pure surface current density as the 
stable equilibrium. 

A. 2 Explicit solutions with minimum magnetic energy 

To further illustrate our theorem we present here three explicit solutions 
for current distributions and magnetic fields that minimize the magnetic 
energy. We do not repeat the solution for a torus since it is a bit lengthy 
and has been published several times already, probably first by Fock [53] . 
but, independently, several times since then, see e.g. [Ml E3 ESI EH ESI EH]- 
Karlsson [12], however, was probably the first to notice that the solution 
minimizes magnetic energy for constant flux. Dolecek and de Launay [15] 
verified experimentally that a type I superconducting torus behaves exactly 
as the corresponding classical perfectly diamagnetic system for field strength 
below the critical field. Here we treat three cases all involving a constant 
external field. For a cylinder, perpendicular to the field, and for a sphere, 
analytical solutions are found. Finally, for a cube with a space diagonal 
parallel to the field, we present a numerical solution. 

A. 2.1 Cylinder in external perpendicular magnetic field 

Consider an infinite cylindrical ideal conductor with radius i? in a external 
constant perpendicular magnetic field. To get the vector potential one must 
solve the following differential equation, 

V X S = V X (V X A) = 0. (61) 

To solve this one should look for the symmetries of the system. We assume 
that the external constant magnetic field points in the y-direction and that 
the cylinder axis coincides with the z-axis. There will then be no dependence 
on the 2;-coordinate so the magnetic field is. 

Moreover, due to the symmetry of the system, the z-component of the mag- 
netic field must be zero, 

|;(P^.)-^-o. (63) 
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These assumptions and constraints transform eq. (IM]) into, 

VxR--(^— ^^ ^-^ l^V-0 
\p^ (9^9^ dp^ p dp j 

which simply is Laplace equation in cylindrical coordinates {p,ip,z). 

Before writing down the general solution, let us consider the boundary 
conditions. As p — )■ cxd, the magnetic field must approach the external one: 
-Bo = B^y = Bolif cos if + psimp). Furthermore, since the magnetic field 
is zero inside the perfect conductor, one concludes from eq. fp21) that the 
vector potential vector must be constant inside the cylinder. Therefore, the 
solution is, 

Az = const. + Bq p cos if, (65) 

\ P J 

for p > R. The magnetic field outside the cylinder becomes, 

B = pBoil ^ I sin v9 + (^ 5o 1 1 + ^- 1 cos v?, (66) 

which implies that, 

B{p = R) = 2Bo cos (p (p (67) 

The magnetic field on the cylinder's surface determines the surface current 
according to eq. (1261) . so we get, 

k = ——BQCOSif z. (68) 

The total current obtained through integration of the surface current is zero 
as expected, otherwise the energy would diverge. A more detailed analysis 
on this problem can be found in |40j . 



A. 2. 2 Superconducting sphere in constant magnetic field 

Similar calculations can be performed for a superconducting sphere with 
radius R in a. constant external magnetic field pointing in the direction of 
the 2;-axis. As for the cylinder case, eq. ( 16T]) is considerately simplified using 
the symmetries of the system. Since the external constant magnetic field 
points along the 2;-axis, there won't be any dependence on the if coordinate 
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and the magnetic field along this coordinate must be zero. Therefore, the 
magnetic field simplifies to, 

B = ^^{A^sme)f--^{A^r)6, (69) 

r sm 6 oB r or 

where we use spherical coordinates r, 9, ip. Again, using these assumptions 
and constraints eq. fl6T|) becomes. 



V X B = -- 
r 



dr [dr ^^ ^V ^ rde Isin^S^ 



(p = (70) 



Since the magnetic field is zero inside the sphere, eq. ( |69|) implies that the 
vector potential has the form, 

A^{r <R) = -^ (71) 

r sm b' 

where C is a constant. To prevent the vector potential from diverging at 
r = and ^ = 0, the constant C must be zero. Furthermore, as r — )■ oo, the 
magnetic field must go to the external field, Bq = BqZ = Bq^v cos 6' — ^ sin 9). 
Therefore, the solution of eq. (17U]) for this case is, 

A^{r>R) = ^(r-^]sme, (72) 

which leads to the following magnetic field outside the sphere, 

f R^\ - f IR^\ 

B = fBo[l--^j cos9-OBo [l + --^j sin^. (73) 

The magnetic field at the sphere surface is thus, 

B = --Bosmee. (74) 

One notes that this is the same field as that of section 13.21 at the surface of 
the inner sphere when energy is minimized. 

Using eq. (I26l) . the surface current density becomes, 

3c 

k = -—Bosm9(p. (75) 

ovr 

Unlike the infinite cylinder in a perpendicular external field, the sphere must 
have a total non-zero electric current, / = ^RBq, to keep the magnetic field 
from entering. A similar approach to this problem can be found in 
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A. 2. 3 Perfectly conducting cube in constant magnetic field 

Starting from the magnetic energy functional we tiave made finite element 
calculations of the current and magnetic field produced when a perfectly 
conducting cube is placed in a (previously) constant field. The constant 
external field is produced by given currents on the surface of a sphere that 
encloses the cube, as in Subsec. 13.11 The cube is placed inside this sphere 
with one of its space diagonals parallel to the external field. The calculations 
verify the results of the variational principle: current density and magnetic 
field become zero inside the cube and the surface current density adjusts to 
achieve this. In Figures and [7] we show results for the current distribution 
of the surface of the cube as seen first from the top (the direction of the 
external field) and then from the side. One observes that the current density 
concentrates along the edges that form a closed path round the cube. 

We also calculated the magnetic fiux through the circular surface enclosed 
by the equator of the enclosing sphere, which we take to have radius R = 1 
and to produce the constant magnetic field B^. = 1 in the interior, when 
empty. For this case the fiux becomes, 

$0 = B^ttR^ = tx^ 3.1416 (76) 

When an ideally conducting sphere of volume V = 1, i.e. radius a = y J;, is 
placed inside (as in Section [3]) the fiux is reduced to. 



<^sp = ^ (l - 1^) ~ 2.3916 (77) 

When the sphere is replaced by a cube of volume V = 1 and a space diagonal 
parallel to the external fields we find the magnetic fiux, 

$cu ~ 2.2733 (78) 

numerically. One notes that such a cube excludes more fiux than a sphere of 
the same volume. 
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Figure 6: Result of finite element calculation of the current distribution on an 
ideally conducting cube in a constant magnetic field. The external field is along a 
space diagonal of the cube and points up from the figure. The current is on the 
surface and is indicated by arrowheads. 
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Figure 7: Same as Figure [6] but with a sideways view of the cube. The current is 
seen to concentrate along the edges that form a closed path round the cube. The 
induced surface current expels the external field from the interior of the cube. 
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